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A CLASS OF MATRIX-VALUED SCHRODINGER OPERATORS 
WITH PRESCRIBED FINITE-BAND SPECTRA 

FRITZ GESZTESY AND LEV A. SAKHNOVICH 


Abstract. We construct a class of matrix-valued Schrodinger operators with 
prescribed finite-band spectra of maximum spectral multiplicity. The corre¬ 
sponding matrix potentials are shown to be stationary solutions of the KdV 
hierarchy. The methods employed in this paper rely on matrix-valued Herglotz 
functions, Weyl-Titchmarsh theory, pencils of matrices, and basic inverse spec¬ 
tral theory for matrix-valued Schrodinger operators. 


1. Introduction 


While basic aspects of inverse spectral theory for matrix-valued Schrodinger op¬ 
erators were established some time ago, finer properties such as isospectral sets 
(manifolds) of potentials, for instance, in the periodic or algebro-geometric finite- 
band cases, are still in their infancy. This paper intends to make a modest con¬ 
tribution to this circle of ideas. More precisely, given a closed set E C 1 of the 
type 


where 


n— 1 


E — | [J [E 2 j,E2j+l]j U [E2n,00), 


( 1 . 1 ) 


{E(} 0 <e< 2 n Cl, n e N, with E t < E e+1 , 0 < i < 2n - 1, (1.2) 

we construct m x m matrix-valued Schrodinger operators Hy = —d 2 /dx 2 I m + Qy 
in L 2 (R) mxm ( I m the identity matrix in C mxm , meN) with prescribed spectrum 

E, 


spec(Rs) = £, (1.3) 

of uniform spectral multiplicity 2m. The constructed matrix potentials Qy will 
turn out to be reflectionless in the sense discussed in [23], [40], and [56], that is, the 
half-line Weyl-Titchmarsh matrices A i±,s(z,x) associated with Hy, the half-lines 
[a;, ±oo), and a Dirichlet boundary condition at x £ R, satisfy 

lim A4 +,e(A + ie,x) = lim A4-.s(A — is, x), (1.4) 

efO e|0 

n —1 

A E [J (E 2 j, E 2j+1 ) U ( E 2n , oo), x E K. 

3=0 


Date : February 1, 2008. 

Key words and phrases. Matrix-valued Schrodinger operators, finite-band spectra, Weyl- 
Titchmarsh matrices. 


1 



2 


F. GESZTESY AND L. A. SAKHNOVICH 


Especially, A^+ i e(-, x) is the analytic continuation of x) through the set E, 

and vice versa. In other words, A4 +,e(-,£) and A4 _,e(-, a;) are the two branches of 
an analytic matrix-valued function A4s(-,x) on the two-sheeted Riemann surface 
of (n?=o( 2 — E i)) 1/2 - These facts imply a purely absolutely continuous spectrum 
E of the associated Schrodinger operator Hy of uniform (maximal) multiplicity 2m. 

Before we turn to a brief description of the contents of each section, it seems 
appropriate to mention some of the pertinent results and especially, the most recent 
activities in connection with (inverse) spectral theory of matrix-valued Schrodinger 
operators. The basic Weyl-Titchmarsh theory of singular Hamiltonian systems and 
their basic spectral theory were developed by Hinton and Shaw, Kogan and Rofe- 
Beketov, Orlov, and others (see, e.g., [43], [45]—[49], [54], [55], [57], [58], [81], [91], 
[92], [94], [97], [101, Ch. 9], and the references therein). Various aspects of direct 
spectral theory, including investigations of the nature of the spectrum involved, 
(regularized) trace formulas, uniqueness theorems, etc., appeared in [10], [12]—[15], 
[20], [23], [24], [35], [37], [39], [60], [82], [89], [90]. General asymptotic expansions 
of Weyl-Titchmarsh matrices as the (complex) spectral parameter tends to infin¬ 
ity under optimal regularity assumptions on the coefficients are of relatively recent 
origin and can be found in [20], [21] (see also [95], [107]). The inverse scattering 
formalism for Schrodinger operators has been studied by a variety of authors and we 
refer, for instance, to [1], [2], [3], [7], [75], [77], [78], [108]. General inverse spectral 
theory, the existence of transformation operators, etc., are discussed in [68], [69], 
[85], [94], [97], [98], [100], [101], and the references therein. Inverse monodromy 
problems have recently been discussed in [4], [5], [6], [15], [68], [69], [96], [101], 
and the literature cited therein. More specific inverse spectral problems, such as 
compactness of the isospectral set of periodic Schrodinger operators [13], special 
isospectral matrix-valued Schrodinger operators, and Borg-type uniqueness theo¬ 
rems (for periodic coefficients as well as eigenvalue problems on compact intervals) 
were recently studied in [17], [18], [19], [21], [23], [26], [51], [52], [68], [69], [101], 
[102], [103], [104], Moreover, direct spectral theory in the particular case of periodic 
Schrodinger operators (i.e., Floquet theory and the like) has been studied in [12], 
[14], [23], [25], [26], [34], [53], [59], [60], [86], [102], [103], [110]-[112], with many 
more pertinent references to be found therein. Apart from Floquet theoretic ap¬ 
plications in connection with Schrodinger operators already briefly touched upon, 
we also need to mention applications to random Schrodinger operators associated 
with strips as discussed, for instance, in [16], [54], [56], and especially to nonabelian 
completely integrable systems. Since the literature associated with the latter topic 
is of enormous proportions, we can only refer to a few pertinent publications, such 
as, [7], [8], [11], [27], [30], [70], [71], [75], [79], [80], [83], [93], [95]-[97], [99], The 
interested reader will find a wealth of additional material in these references. 

Section 2 summarizes basic results in Weyl- Titchmarsh theory and some ele¬ 
ments of inverse spectral theory for matrix-valued Schrodinger operators. Polyno¬ 
mial pencils of matrices are briefly reviewed in Section 3. In Section 4 we present 
our principal new result, the construction of m x m matrix-valued Schrodinger op¬ 
erators Hy with spectrum E (cf. (1-1)) and uniform maximal spectral multiplicity 
2m. In our final Section 5 we prove that Qs satisfies a stationary KdV equation 
(in fact, we explicitly identify the first equation in the stationary KdV hierarchy 
satisfied by Qy) and derive matrix-valued trace formulas for Qy and higher-order 
KdV invariants. 
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2. Basic Facts on Weyl-Titchmarsh Theory 


In this section we briefly recall basic elements of the Weyl-Titchmarsh theory 
for matrix-valued Schrodinger operators. Throughout this paper all matrices will 
be considered over the field of complex numbers C, and the corresponding linear 
space of k x i matrices will be denoted by C fcx ^, k,£ £ N. Moreover, Tj- denotes the 
identity matrix in C kxk , M* the adjoint (i.e., complex conjugate transpose), M* 
the transpose of a matrix M, diag(mi,... , to*,) £ C fcxfe a diagonal k x k matrix, 
and ACioc(R) denotes the set of locally absolutely continuous functions on R. The 
spectrum, point spectrum (the set of eigenvalues), essential spectrum, absolutely 
continuous spectrum, and singularly continuous spectrum of a self-adjoint linear 
operator T in a separable complex Hilbert space are denoted by spec(T), spec p (T), 
speCes S (T), spec ac (T), spec sc (T), respectively. 

The basic assumption for this section will be the following. 


Hypothesis 2.1. 

(i) Fix m £ N, suppose Q = Q* € L 1 1 oc (R) mxm and introduce the differential 
expression 

d 2 

r = -T m ^ + Q, (2.1) 

(ii) Suppose C is in the limit point case at ±oo. 


Given Hypothesis 2.1 (i) we consider the matrix-valued Schrodinger equation 

—i(>"(z,x) + Q(x)if>(z,x) = zip(z,x) for a.e. xtl, (2.2) 


where z € C plays the role of a spectral parameter and if is assumed to satisfy 


if>{z, •), '•P'iz, •) G AC loc (R) mxm . (2.3) 


Throughout this paper, ^-derivatives are abbreviated by a prime /. 

Let ty(z,x,x o) be a 2m x 2 to normalized fundamental system of solutions of 
(2.2) at some xq £ R which we partition as 


^{z,x,x 0 ) 


fO(z,x,x 0 ) (f)(z,x, Zo)\ 

\0'(z,x,x o) <f)'{z,X,Xo)J 


(2.4) 


Here / denotes d/dx, 9(z,x,x o) and <f>{z,x,x o) are to x to matrices, entire with 
respect to z £ C, and normalized according to 


T(z, a; 0 ,a;o) =F 2m , 


(2.5) 


that is, 

0(z,x o ,x o ) = 4>'(z,x 0 ,x 0 ) =l m , (f>(z,x 0 ,x 0 ) = 9'(z,x 0 ,x 0 ) = 0. (2.6) 

In this context, we briefly recall a set of formulas needed later in Section 4. Intro¬ 
ducing 

J = (l ~o m ) ■ « 2 - 7 ' 

one infers 


^(z,x,x 0 y J^{z,x,x 0 ) = J, (2.8) 

which implies J^{z, x, xq){^{z, x, xq)J)* = T 2m and hence 

ty(z,x,xo)J'S’{z,x,x 0 )* = J ■ 


(2.9) 
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Writing out (2.8) and (2.9) explicitly yields 

9'{z , x, xo)*9(z , a;, xo) — 9(z, x, xo)*9’(z, x, xo ) = 0, (2-10) 

4>{z, x, x 0 )*<f>(z, x , xo) - x , x 0 )*<j>'(z, x , Xo) = 0, (2.11) 

</>'(£, a:, aio)* 0 ( 2 , z, x 0 ) - <f>{z, x, x o )*0'(z, x, x 0 ) = J m , (2-12) 

9{z, x, x 0 )*<j)'(z, x, x 0 ) - 0'(z, x, xo)*<j>{z, x, x 0 ) = T to , (2-13) 

and 

4>{z, x, xo)9(z, x, Xo )* — 0 ( 2 , x, xq)4>{z, x, Xo)* = 0, (2-14) 

<f>'{z, X, Xo )9'(z, x, x 0 )* - 0'(z, x, x 0 )4>{z, x, x 0 )* = 0, (2-15) 

<fi'(z, x, x 0 )9(z, x, x 0 )* - 9'(z, x, x 0 )<t>{z, x, x 0 )* = l m , (2.16) 

9(z, x, x 0 , a)<t>'(z, x, x 0 )* - 4>{z, x, x 0 , a)9'(z, x, x 0 )* = I m . (2.17) 

Next, assuming —00 < a < b < 00 , we consider the spaces 


N(z, ± 00 ) = {(j> £ L 2 {(c , ±oo)) m | —(f>" + Q<fi = z<p a.e. on (c, ± 00 )}, (2-18) 

for some c6l and z £ C. (Here {<f>, iff) o = J2j =1 denotes the standard scalar 
product in C™, abbreviating y € C™ by x = (Xi; • ■ ■ , Xu)*-) Both dimensions of the 
spaces in (2.18), dimc(-/V(z, 00 )) and dimc(./V(z, — 00 )), are constant for z £ C± = 
{£ G C | ±Im(£) > 0} (see, e.g., [55]). One then recalls that C in (2.1) is in the 
limit point case at ±00 whenever 

dime (AT ( 2 , ± 00 )) = m for all 2 € C\R. (2-19) 

Since the potential Qs to be constructed in Section 4 will automatically lead to 
the limit point case at ± 00 , we decided to limit our considerations mainly to this 
situation. In this context we note the well-known fact that if C in (2.1) is in the 
limit point case at ± 00 , then the mxm Weyl-Titchmarsli matrices associated with 
£, the half-lines [x, ± 00 ), and a Dirichlet boundary condition at x, are given by 

M±(z,x 0 ) = y , ± {z,x,xo)'&±(z,x,xo)- 1 \ x=xo , z £ C\R, (2.20) 
where satisfy (£ — zl m )'f>±(z, -,xo) = 0 and 

*±(z, ; x 0 ) £ L 2 ([x 0 , ± 00 ))“. (2.21) 

The actual normalization of \i/±(z, -,Xq) is clearly irrelevant and hence ^±( 2 , -,Xo) 
can be replaced by ^f±(z, -,Xo )C, where C is any nonsingular mxm matrix. 

For later reference we summarize the principal results on A4±(z,Xo) in the next 
theorem. First we recall the following definition. 

Definition 2.2. A map A4: C+ —*• C nxn , n £ N, extended to C_ by A i(z) = 
A i(z)* for all z £ C+, is called an n x n Herglotz matrix 1 if it is analytic on C+ 
and Im(A4(z)) > 0 for all z £ C+. 

Here we denote Im(Af) = (At — M*)/2i and Re(A4) = (At + At*)/2. 

In the scalar context n = 1, the condition Im(At(z)) > 0 in Dehnition 2.2 can 
be replaced by Im(m(z)) > 0 for the corresponding scalar counterpart m{z). 

1 There appears to be considerable confusion in the literature since Nevanlinna, Pick, 
Nevanlirma.-Pick matrix, in addition to Herglotz matrix, are also in use. In part these discrepan¬ 
cies can be traced back to the use of the upper half-plane C+ versus the open unit disk D, and in 
some cases the geographical location of the author in question determines the preferred notation. 
Following a tradition in mathematical physics, we adopt the notion of Herglotz functions in this 
paper. 
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Theorem 2.3 ([21], [40], [45], [46], [49], [56]). Assume Hypothesis 2.1 and suppose 
that z £ C\R, xo £ R. Then 

( i ) ±A4±(z,Xq) is an m x m matrix-valued Herglotz function of maximal rank. In 
particidar, 


lvn(±M.±{z,X(f)) > 0, 2 € C+, 

M±(z,x 0 ) = M±(z,x 0 )*, 
rank(A4±(2, Xo)) = m, 
lim _A4±(A + ie, xo) exists for a. e. A £ R. 

ej.0 


( 2 . 22 ) 

(2.23) 

(2.24) 

(2.25) 


Local singularities of ±A4±(z, xo) and tA4±(z, xq) 1 are necessarily real and at 
most of first order in the sense that 

(2.26) 


=F lim (ie A4±(A + ie, xo)) > 0, A £ R, 

e|0 

± limie.A/(±(A + ie, Xo) -1 > 0, A £ M. 

e|0 


(2.27) 


(ii) ±J\4±(z, xo) admit the representations 

±M±(z,x 0 ) = Re(±A4±(±i,xo)) + f dfi±(A,x 0 ) ((A - z) -1 - A(1 + A 2 ) -1 ), 

Jr 


where 


and 


I (A, Xo) ||c mx ' m (1 A ) < oo 

Jr 


(2.28) 

(2.29) 


^ rn+d 

fi±((A, /a], xo) = limlim — / dvlvn(izM±{v + ie, xo)). (2.30) 

<5X0 eXO 7T J x+S 


(Hi) Define the 2to x m matrices 

^±(z,x,x 0 ) 


'&±(z,x,x o) = 


ip±(z,x,x o) 

9(z,x,x 0 ) <f>(z, x, Xo) 

9'(z, X, X 0 ) <j>'(z, x, x 0 )) \M±(z,x 0 )J ’ 


(2.31) 


then 


/*±oo 

Im(A4±(-2, Xo)) = Im(z) / dx if±(z, x, xo)*t/>±( 2 , x, xo). (2.32) 
J Xo 

(iv) Denote by C e C C+ t/ie open sector with vertex at zero, symmetry axis along 
the positive imaginary axis, and opening angle e, with 0 < e < tt/ 2. Then 


M±(z,x 0 ) = ±z 1/2 l m + o(l). 

\z \—»oo 

zec £ 


(2.33) 


Necessary and sufficient conditions for M±(z, Xo) to be the half-line toxto Weyl 
Titchmarsh matrix associated with a Schrodinger operator on [xo,±oo) in terms 
of the corresponding measures fl±(-,Xo) in the Herglotz representation (2.28) of 
M±(z, Xo) can be derived using the matrix-valued extension of the classical inverse 
spectral theory approach due to Gelfand and Levitan [33], as worked out by Rofe- 
Beketov [85]. The following result describes sufficient conditions for a monotonically 
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nondecreasing matrix function to be the matrix spectral function of a half-line 
Schrodinger operator. It extends well-known results in the scalar case to = 1 (cf. 
[65, Sects. 2.5, 2.9], [ 66 ], [76, Sect. 26.5], [106]). 

Theorem 2.4 ([85]). Let be a monotonically nondecreasing to x m matrix¬ 
valued function on R satisfying the following two conditions. 

(■ i) Whenever f G C([xo, oo)) mxl with compact support contained in [xo,oo) and 


where 


(ii) Define 


[ F(A)*dQ+(X) F(\) = 0, then f = 0 a.e., 

Js. 


.. f R , sin(A 1 / 2 (a; - xo)) , . 

F(A) = hm / dx - 7777 ;-Ag 


R t oo t 


A 1 / 2 


n fA) = K( A )-Jf A3/2 - A ^° 

+ ( } |fi+(A), A < 0 


and assume the limit 


lim 

R'loo 


dfl+{ A) 


sin(A 1 / 2 (a; — xo)) 
AV2 


$(a;) 


(2.34) 


(2.35) 


(2.36) 


(2.37) 


exists and 4> G L°°(\xq, R\) mxm for all R > xo■ Moreover, suppose that for some 
r G No, 4>( r+1 ) G L\[xq, R]) mxm for all R > xq, and $(xo) = 0. 

Then 17+ is the matrix spectral function of a self-adjoint Schrodinger operator H + 
in L 2 ([x o,oo)) m associated with the m x to matrix-valued differential expression 
£+ = —d 2 /dx 2 T m + Q, x > x o, with a Dirichlet boundary condition at xo, a self- 
adjoint boundary condition at oo (if necessary), and a self-adjoint potential matrix 
Q with Q^ G L 1 ([xo,R\) rnxrn for all R > x 0 . 


Next, assuming Hypothesis 2.1, we introduce the self-adioint Schrodinger oper¬ 
ator H in L 2 (R) m by 

d 2 

H = (2.38) 

dom(ir) = {g G L 2 (R) m | 5 , 5 ' G AC loc (R) m ; (- g" + Qg) G L 2 (R) m }. 


The resolvent of H then reads 

((H - 0 )" 1 /)( x) = f dx' Q(z, x, x')f(x'), z G C\R, / G L 2 ( R) m , (2.39) 
it 


with the Green’s matrix Q(z, x, x') of H given by 

Q(z,x,x') = ip T (z,x,x 0 )[M-(z,xo) - M+(z, x 0 )] _ V±(^, x',x 0 )*, 

i^i'.zG C\R. (2.40) 

Introducing 


Jf±(z,x o) = M-(z,x 0 ) ± M + (z,x 0 ), z G C\R, 


(2.41) 



MATRIX SCHRODINGER OPERATORS WITH FINITE-BAND SPECTRA 


7 


the 2 to x 2m Weyl-Titchmarsli function M(z, xq) associated with H on R is then 
given by 


M(z,x 0 ) = o )) jy=12 


fM±(z, x 0 )N-(z,xo) 1 M t {z,xq) 
\ M+{z,x 0 )N-{z,x 0 )- 1 /2 


(2.42) 



The basic results on M(z, xq) then read as follows. 


Theorem 2.5 ([40], [45], [46], [49], [56]). Assume Hypothesis 2.1 and suppose that 
z £ C\R, xq £ R. Then, 

(i) M(z,x o) is a matrix-valued Herglotz function of rank 2m with representation 


M(z, xo) = Re(M(i,xo)) + ( dCl(X, x 0 ) ((A - z) 1 - A(1 + A 2 ) x ), (2.43) 

Jr 

j ||dfi(A,£o)||c 2mx2m (1 + A 2 ) -1 < oo (2.44) 

Jr 

^ rU +<5 

f2((A, p], xq) = limlim — / dulm(J\4(u + ie, xo))- (2.45) 

<5t0 eiO 7T J x+S 

( ii ) 2 £ C\spec (H) if and only if M(z, xq) is holomorphic near z. 


where 


and 


Here spec(T) denotes the spectrum of T. Later on we will denote by spec ac (T) 
the absolutely continuous spectrum of T. 

Finally, we state the following characterization of M(z,Xo) to be used later on. 
In the scalar context m = 1 this has been used by Rofe-Beketov [871, [881 (see also 
[65, Sect. 7.3]). 

Theorem 2.6 ([87], [ 88 ]). Assume Hypothesis 2.1, suppose that z £ C\R, Xo € R, 
and let f,r£No. Then the following assertions are equivalent. 

(i) Ai(z,Xo) is the 2mx2m Weyl-Titchmarsh matrix associated with a Schrodinger 
operator H in L 2 (R) m of the type (2.38) with an to x to matrix-valued potential 
Q £ Ll oc (R) and Q £ C e ((—oo,x 0 )) and Q £ C r ((x o,oo)). 

(ii) A4(z,xq) is of the type (2.42) with A / t±(z,xo) being half-line m x m Weyl- 
Titchmarsh matrices on [xo, ±oo) corresponding to a Dirichlet boundary condition 
at xq and a self-adjoint boundary condition at — oo and/or oo (if any) which are 
associated with an m x m matrix-valued potential Q satisfying Q £ C e ((—oo,Xo)) 
and Q £ C r ((X q,oo)), respectively. 

If (i) or (ii) holds, then the 2m x 2m matrix-valued spectral measure fl(-,x o) asso¬ 
ciated with A4 (z,Xq) is determined by (2.42) and (2.45). 


Next, we consider variations of the reference point xq £ R. In analogy to (2.20), 
we note that in the case where the Schrodinger differential expression £ is in the 
limit point case at ±oo, 

M±(z,x) = ^' ± x (z,x, xo)'&±(z,x, xo) -1 , z £ C\R, (2.46) 

represents the corresponding half-line Weyl-Titchmarsh matrix on [x, ±oo), iel, 
with ^±(z, •, xo) defined in (2.31). Again the actual normalization of \H± is, of 
course, irrelevant. Since satisfies the second-order linear m x m matrix-valued 
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differential equation (2.2), Al± in (2.46) satisfies the matrix-valued Riccati-type 
equation (independently of any limit point assumptions at ±oo) 

M'±(z, x) + M±(z, x) 2 — Q(x) — zl m , x G R, z G C\M. (2.47) 

The asymptotic high-energy behavior of A i±(z,x) as \z\ —> oo has recently 
been determined in [21] under minimal smoothness conditions on Q and without 
assuming that C is in the limit point case at ±oo. Here we recall just a special case 
of the asymptotic expansion proved in [ 21 ] which is most suited for our discussion 
in Section 4. We denote by C e C C+ the open sector with vertex at zero, symmetry 
axis along the positive imaginary axis, and opening angle e, with 0 < e < 7 t/2 . 

Theorem 2.7 ([21]). Fix Xq € R and let x > xq. In addition to Hypothesis 2.1 
suppose that Q G C°°([a;o, ±oo)) mxm and that C is in the limit point case at ±oo. 
Let A4±(z,x), x > xo, be defined as in (2.46). Then, as |z| —> oo in C e , A i±(z,x) 
has an asymptotic expansion of the form (Im(^ 1 / 2 ) > 0 , z £ C+) 

N 

M±(z,x) = ±il m z 1 ^ 2 + M±,k(x)z~ k ^ 2 + o(|z|~ jY/2 ), Ng N. (2.48) 

\z \—►OO ‘ ^ 

1 1 7 1 


The expansion (2.48) is uniform with respect to arg( 2 ;) for \z\ —> oo in C £ and 
uniform in x as long as x varies in compact subsets of [xo,oo). The expansion 
coefficients A4±,fc(x) can be recursively computed from 

•Adi,i(x) =;*||Q(x), A1±, 2 (x) = ^Q'(x), 

. , k -1 x (2.49) 

Al±,fc + i(x) = M± tk (x) + ^ M± 1 e(x)M±,k-e(x) j , k > 2 . 

^ 7-1 ' 


The asymptotic expansion (2.48) can be differentiated to any order with respect to 
x. 


Remark 2.8. 

(i) Due to the recursion relation (2.49), the coefficients M±,k are universal poly¬ 
nomials in Q and its x-derivatives ( i.edifferential polynomials in Q). That the 
asymptotic expansion (2.48) can be differentiated to arbitrary order in x follows 
from repeated use of the Riccati-type equation (2.47). 

(ii) In the case where Q and its x-derivatives are in L 1 (K) mxm , or in the case 
where Q is periodic and hence Floquet theory applies, the proof of the existence of 
an asymptotic expansion of the type (2.48) follows in a routine manner by iterating 
appropriate Volterra-type integral equations. The general case, however, is intricate 
as is evident from the treatment in [ 21 ]. 

Finally, in addition to (2.40), one infers for the 2m x 2m Weyl-Titchmarsh 
function M(z,x) associated with H on R in connection with arbitrary half-lines 



MATRIX SCHRODINGER OPERATORS WITH FINITE-BAND SPECTRA 9 

[x, ±oo), X € R, 

M(z,x) = (Mj^{z,x)) jjl=12 , zeC\K, (2.50) 

x) = M±(z, x)[M-{z, x) - M+(z, x)]~ 1 M T {z, x) 

= ^)' + {z,x,xq)[M-{z,x 0 ) - M+(z,x,x 0 y, (2.51) 

Mi,2(z, x) = 2 ~ 1 [M~(z, x) - M+(z, x)]~ 1 [M-{z, x) + M+(z, x)} 

= ip + {z,x,x 0 )[M-{z,x 0 ) - M+(z, xo)] -1 ^ (z,x,xo)*, (2.52) 

M2,i(z, x) = 2 ~ 1 [M-(z, x) + M + (z, x)][M-(z, x) - M+(z, a ;)] -1 

= ip' + (z,x,x 0 )[M-(z,x 0 ) - M + {z,x 0 )}~ 1 ip-(z,x,x 0 y, (2.53) 

M2,2(z, x) = [M~{z, x) - M+(z, a;)]^ 1 

= i{>+(z,x,xo)[M-(z,xo) - M + (z,xo)]~ 1 ip-(z,x,x 0 y. (2.54) 

Introducing the convenient abbreviation, 

"(>.*) = (_%% -££?)■ *■» 

one then verifies from (2.50)-(2.55) and from M(z,x)* = M(z,x), M±(z,x)* = 
M±(z,x) that 


g(z,x)* =g{z,x), g 2 (z,x)* = Qi(z,x), f)(z,x)* = t){z,xfi (2.56) 

g(z,x)gi(z,x) = g 2 (z,x)g(z,x), (2.57) 

f )(z, x)g 2 (z, x) = 0 i( 2 , ar)h (z, x), (2.58) 

&(z, x) = [M-(z, x) - M+(z , a;)] -1 , (2.59) 

9(z, x)f)(z, x) - g 2 (z, xf = -(l/4)T m , (2.60) 

f )(z,x)g(z,x) - gi(z,x) 2 = -(l/4)T m , (2.61) 

M±(z, x) = =F(l/2)g (z, a ;) -1 - g(z, x)~ 1 g 2 (z, x) (2.62) 

=F (l/ 2 ) 0 ( 2 , a ;) -1 - 01 ( 2 , x)g(z, a;) -1 . (2.63) 


Moreover, the Riccati-type equations (2.47) imply the following results needed in 
Section 4. 


Lemma 2.9. Let z £ C\R and define M± by (2.46) so that M± satisfy the Riccati- 
type equation (2.47). Then, for a.e. x £ R, 


0 ' =-( 01 + 02 ), (2-64) 

g'i = -{Q-zl m )g-\) (2.65) 

= (-0" + 0Q-<20)/2, (2-66) 

02 = ~&(Q ~ zTm) - f) (2.67) 

= (-0" + Q0 - 0(2)/2, (2.68) 

f)' = - 0 i (Q - zT m ) - (Q - zl m )g 2 , (2.69) 

f) = [0" - 0(Q - zlm) - (Q - ^m)0]/2 (2.70) 
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if Q € Ll oc (JSL) mxm , and 

01 = -2 (Q - zl m ) g ' - Q'd + 01Q - 201, (2.71) 

02 = -20'(Q - zlm) - 02' + 202 - 02 2 (2.72) 

if in addition Q! £ L 1 1 oc (IR) mxm . 

Proof. (2.47) rewritten in terms of 0, 0i, 02 yields 

± ( 1 / 2 ) 0 _ 1 0 , 0 1 + 0 _ 1 0 / 0 _1 02 - 0 _1 02 + ( l / 4 ) 0^ 2 + 0 ~ 1 020~ 1 02 

± (l/2)0- 2 0 2 ± (l/2)0- 1 0 2 0- 1 = Q - zl m . (2.73) 

Taking the difference of the two equations in (2.73) yields (2.64). Adding the 
two equations in (2.73) and using (2.57), (2.60), (2.61), and (2.64) yields (2.65) 

and (2.67). Combining (2.65), (2.67), and (2.64) implies (2.70). Inserting (2.70) 

into (2.65) and (2.67) yields (2.66) and (2.68). (2.69) follows from differentiating 
f) = 010' 1 — (1/4)0 _1 , inserting 0' and 0) from (2.64) and (2.65), and making re¬ 
peated use of the identities (2.57), (2.61). Finally, (2.71) (resp. (2.72)) follows from 
differentiating (2.65) (resp. (2.67)) inserting (2.69) for (/. Alternatively, (2.67)- 
(2.72) follow directly from (2.65)-(2.71) using (2.56). □ 

3. Polynomial Pencils of Matrices in a Nutshell 

Since self-adjoint polynomial pencils of matrices play a role in our principal 
section 4, we briefly review some of the corresponding definitions and basic results, 
mainly following the monograph of Markus [72] and papers by Markus and Matsaev 
[73], [74], While all results below are discussed for operator pencils by Markus and 
Matsaev, we will only quote them in the matrix context, for simplicity. 

Given m € N, we denote by 

n 

A(z) = J2 A kz k , A k e C mxm , 1 < k < n, z £ C, (3.1) 

k -0 

a polynomial pencil of m x m matrices (in short, a pencil) in the following. A is 
called of degree n € No if A n ^ 0 and monic if A n = T m . 

Definition 3.1. Let A be a pencil of the type (3.1). 

(f) The spectrum of A, denoted by spec(A) , is defined by 

spec(A) = {z £ C | A{z) is not invertible}. (3.2) 

zo € C is called an eigenvalue of A if A(zo)fo = 0 has a solution fo € C m \{0}. 

( ii) A monic pencil C is called a (right) divisor of A if 

A(z) = B(z)C(z), zGC (3.3) 

for some pencil B. If in addition spec(£>) D spec(C) = 0, then C is called a (right) 
spectral divisor of A. 

(Hi) Z £ C mxm i s called a (matrix) root of the pencil A if 

A(Z) = 0, (3.4) 

where A(Z) is defined as 

n 

A(Z) = Y, A k zk - 

k—0 


(3.5) 
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Z G C mxm is called a ( matrix) spectral root of the pencil A if (zT m — Z) is a 
spectral divisor of A. 

(iv) The pencil A is called self-adjoint if Ak = Af, for all 1 < k < n ( i.e., A(z)* = 
A(z) for all zeC). 

(v) A self-adjoint pencil A is called weakly hyperbolic if A n > 0 and for all f G 
C m \{0}, the roots of the polynomial (f,A(-)f)c m are real. If in addition all these 
zeros are distinct, the pencil A is called hyperbolic. 

(vi) Let A be a weakly hyperbolic pencil and denote by {Pj(A, f)}\<j< n , 

Pj (A,f)<p j+1 (AJ), l<j <n-l, /GC m \{0}, (3.6) 

the roots of the polynomial (f, A(-)f)c m ordered in magnitude. The range of the 
roots Pj(A , /), / G C m \{0} is denoted by A.,-(-4) and called the jth root zone of A. 

(vii) A hyperbolic pencil A is called strongly hyperbolic if A j (A) and Afc(-4) are 
mutually disjoint for j ^ k, 1 < j, k < n. 


Moreover, the Vandermonde matrix corresponding to a collection {Z \,..., Z n } C 
C mxm is defined by 


V(Z 1 ,...,Z n ) 


T 

m 

T 

- L 'm • • 

Zm \ 

Zi 


z n 

jn—1 

-'l 

gm — 1 
^ 2 

. Z n ~ x ) 

n / 


(3.7) 


Theorem 3.2 (Markus [72], Sect. 29). Let A be a pencil of the type (3.1). 

(z) ( zl m — Z) is a divisor of A if and only if A(Z) = 0. ( This justifies the notation 
introduced in the last part of Definition 3.1 ( ii ).) 

(ii) Let A be a monic pencil of degree n and Z\,... ,Z n spectral roots of A. Then 
the following assertions are equivalent: 

(a) V(Zi ,..., Z n ) is invertible. 

(/3) spec (Zj) n spec(Zfc) = 0, j ^ k, 1 < j, k < n. 

(7) spec(M) = Uj=i spec(-Zj). 

(Hi) If A is a self-adjoint pencil, then spec(M) is symmetric with respect to R. 


Theorem 3.3 (Markus [72], Sect. 31). Let A be a pencil of the type (3.1). 

(i) If A is a weakly hyperbolic pencil, then spec(M) C R. 

(ii) The root zones Aj(M), 1 < j < n, of a weakly hyperbolic pencil A are intervals 
(possibly degenerating to a point ). 

(Hi) If A is a weakly hyperbolic pencil of order n with root zones {A J (M)}i< :/ < n 
and Ak(A) = {A k} for some k G {1,..., n}, then A(z) = (z — A k)B(z), where B is 
a weakly hyperbolic pencil of order n — 1, with root zones 

Ai(M),..., Afc_i(M), Afc + i(M),..., A„(M). (3.8) 

(iv) If A is a weakly hyperbolic pencil then Aj(M)nAj + i(M), 1 < j < n— 1, consists 
of at most one point. 

(v) If A is a hyperbolic pencil, then Aj(M)nAfc(M) = 0 , j k, 1 < j, k < n. Thus, 
a hyperbolic pencil is strongly hyperbolic if and only if 

Aj(A) n A J+ i(M) = 0, l<j<n-l. (3.9) 

(vi) Suppose A is a self-adjoint pencil of degree n, A n > 0, and -4(A) ^ 0 for all 
A G R. Then A is a weakly hyperbolic pencil if and only if there exist numbers 
7i < 72 < • • • < 7n-i such that (—iyA^j) >0, 1 < j < n — 1. 
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(vii) Suppose A is a self-adjoint pencil of degree n with A n > 0. Then A is a 
strongly hyperbolic pencil if and only if there exist numbers 71 < 72 < • • • < Jn-i 
such that (— iyA(pij) >0, 1 < j < n — 1 . 

Theorem 3.4 (Markus [72], Sect. 31). Let A be a pencil of the type (3.1). 

(■ i ) Suppose A is a weakly hyperbolic pencil and A J 0 _ 1 (M) D A J 0 (7l) = Aj 0 (A) H 
A j 0 + i(7l) = 0 for some jo £ {1, ... ,n). Then A has a spectral root Zj 0 such that 
spec(A) n A j 0 (A) = spec (Zj 0 ) and ^jo is similar to a self-adjoint matrix. 

(■ ii ) A strongly hyperbolic pencil A has n spectral roots {Zj}\<j< n such that spec(2j) 
= spec(^4) fl Aj(A) and each Zj, 1 < j < n, is similar to a self-adjoint matrix. 

Theorem 3.5 (Markus [72], Sect. 31, Markus and Matsaev [73], [74]). Let A be a 
pencil of the type (3.1). 

(i) A weakly hyperbolic monic pencil A of degree n is decomposable as 

A(z) = (zlm - y n ){zlm ~ y n - 1 ) • • • {zlm ~ ^l), (3.10) 

with spec((Vj) C Aj(A), 1 < j < n. 

(ii) Let A be a strongly hyperbolic monic pencil of degree n. Then A is decomposable 
as 


A(z) = (zT m - y n )(zl m - y n—l) ■ ■ ■ (zlm ~ (Vl ), (3.11) 


with spec(3 / ]j) C Aj(A), 1 < j < n. Moreover, each (Vj is similar to a spectral root 
Zj of A and hence, 

spec((Vj) = spec(Zj) = spec(^l) H Aj(^l), 1 < j < n. (3-12) 


4. A Class of Matrix-Valued Schrodinger Operators 
with Prescribed Finite-Band Spectra 

This section is devoted to the construction of a class of matrix-valued Schrodinger 
operators with a prescribed finite-band spectrum of uniform maximum multiplicity, 
the principal result of this paper. 

To begin our analysis we start with a useful result on (scalar) Herglotz functions. 
Even though the result is probably well-known to experts, we provide an elementary 
proof for completeness. 

Let 


{Ee}o<e< 2 n CM, n e N, with Et < Ei+ 1, 0 < £ < 2n — 1, 
and introduce the polynomial 

2 n 

Rln+\(z) = - El), Z G C. 

1=0 


Moreover, we define the square root of i?2n+i by 

R-2n+ i (^) 1 ^ 2 = hmiWi(A + is ) 1 ^ 2 , A 6 R, 

e|0 


(4.1) 


(4.2) 


(4.3) 
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and 


for A € (— 00 , E 0 ), 
for A G (E 2 j-i,E 2 j), j = 1,..., n, 
for A G (E 2j ,E 2j+ 1 ), j = 0,..., n - 1, 
for A G (E 2n , oo), 

A G R (4.4) 

1 /2 

and analytically continue R 2n+1 from R to all of C\£, where £ is defined by 

n -1 'j 

[J [E 2 j , E 2 j _)_i] > U \E 2n , oo). (4-5) 

i=o J 

In this context we also mention the useful formula 

= -R 2n +i{z) 1/2 , z G c+. (4.6) 



#2n+l(A) 1/2 — |i?2n+l(A) 1/2 | 


Theorem 4.1. Let 2 G C\£ and n G N. Define i? 2 n+i as (4.1)-(4.4) followed 
by an analytic continuation to C\£. Moreover let F n and H n+ 1 6e two monic 
polynomials of degree n and n + 1, respectively. Then 


iF n (z) 


(4.7) 


is a Herglotz function if and only if all zeros of F n are real and there is precisely 
one zero in each of the intervals [E 2 j_\,E 2 j\, 1 < j < n. Moreover, if iR 2n + l F n is 
a Herglotz function, then it can be represented in the form 


iFn(z) = 1 f F n (X)dX 1 
R 2 n+i(zy/ 2 7r7 E i?2 n+1 (A)V2 A-*’ X 


(4.8) 


Similarly, 


iH n+1 (z) 

R 2 n +i(z ) 1/2 


(4.9) 


is a Herglotz function if and only if all zeros of H n+ \ are real and there is precisely 
one zero in each of the intervals (—oo,i?o] and [E 2 j_\,E 2 j], 1 < j < n. Moreover, 
if iR 2n \_ 1 H n+ i is a Herglotz function, then it can be represented in the form 


iH n+1 (z ) / iH n+1 (i) A 

R 2n+ i(z) 1 / 2 ~ 6 \ Rln+l (i) 1/2 / 

1 f H n+ i(X)dX (l X \ 

' n ]■£ R 2 n+ 1 (X) 1 ^ 2 \ A — Z 1 + X 2 )' 


z G C\£. 


(4.10) 


Proof. We start with the case of F n (z) / R 2n+ i(z) 1 ^ 2 in (4.7). Consider a closed 
counterclockwise oriented contour which consists of the semicircle C e = {z G 
C | z = E 0 +e exp(ia), —7 t/2 < a < 7t/ 2} centered at E 0 , the straight line L + = {z G 
C + | z = x + is, E 0 < x < R}, the following part of the circle of radius ( R 2 + e 2 ) 1 / 2 
centered at E 0 , Cr = {z G C | z = E 0 + (R 2 + e 2 ) 1 / 2 exp (i/3), arctan(e/i?) < (3 < 
2n — arctan(e/iZ)}, and the straight line L_ = {zG C_ | z = x — ie, E 0 < x < R}. 
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Then, for £ > 0 small enough and R > 0 sufficiently large, one infers 

iK(z) = J_ / 1 iF n ( C) 

iW(*) 1/2 2mf rRe C-zR 2n+1 (<:Y /2 ^ 

= if — _ F ™( x ) dx u n \ 

eiO.flToo 7T A - 2: i?2n+l(A) 1 /2 

Here we used (4.4) to compute the contributions of the contour integral along [E 0 , i?] 
in the limit e J, 0 and note that the integral over Cr tends to zero 

Next, utilizing the fact that F n is monic and using (4.4) again, 
F n (A)dA/i? 2 n+i(A) 1 / 2 represents a positive measure supported on 
F n has precisely one zero in each of the intervals [E 2 j- 1 , E 2 j\, 1 < 
words, 


as R t 0 since 

(4.12) 

one infers that 
£ if and only if 
j < n. In other 


F n ( A) 

iW(A)V2 


> 0 on S 


(4.13) 


if and only if F n has precisely one zero in each of the intervals [E 2 j-i, E 2 j], 1 < 
j < n. The Herglotz representation theorem, Theorem 2.3, then finishes the proof 
of (4.8). The proof of (4.10) follows along similar lines taking into account the 
additional residues at ±z inside Tr, e which are responsible for the real part on the 
right-hand side of (4.10). □ 


Theorem 4.1 can be improved by invoking ideas developed in the Appendix of 
[61] (cf. also [105]). We will pursue this further in [9]. 

Corollary 4.2. Let z G C\£ and m,n € N. Define as (4.1)-(4.4) fol¬ 

lowed by an analytic continuation to C\£. Moreover let T n and H n +1 be two monic 
m x to matrix pencils of degree n and n + 1, respectively. Then ( i/2)R 2n l l T n is a 
Herglotz matrix if and only if the root zones A j(tF n ) of T n satisfy 

A j{T n ) C [E 2 j-i,E 2 j\, 1 < j < n. (4.14) 

_1 /Q 

Analogously, (i/2)R 2n ^_ 1 7i n +i is a Herglotz matrix if and only if the root zones 
Aj(H n+ i) ofHn+i satisfy 

Ao(W n +i) C (—oo, Ao], Aj(H n+ i) C [E 2 j_i, E 2 j], 1 < j < n. (4.15) 

If (4.14) ( resp (4.15)) holds, then T n (resp ., TL n +i ) is a strongly hyperbolic pencil. 

Proof. We recall that ( i/2)R 2n + 1 tF n is an in x m Herglotz matrix if and only if 

(/, (i/2)Rf R f 2 1 iF n f)c m is a Herglotz function for all / € C m \{0}. Thus, it suffices 
to apply Theorem 4.1, identifying F n and {f,F n f)c m , to arrive at (4.14). The same 
argument applied to H n +1 yields (4.15). □ 

Next, we define the following 2 to x 2to matrix A4s(z, xo) which will turn out to 
be the underlying Weyl-Titchmarsh matrix associated with a class of to x to matrix¬ 
valued Schrodinger operators with prescribed finite-band spectra. We introduce, 
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for fixed Xq G R, 

Ms(z,x o) = (Ats,p,g(^,a;o)) 1 < P)g < 2 (4.16) 

i ( U n +l,Y.{z,X 0 ) -G2,n-l,T.(z,Xo)\ 

2R 2n +l(z) 1/2 V - ^l.n-l.s(^)*o) Fn,-z{z, Xo) ) ’ ' 

Here i? 2 n+i(^) 1 / 2 is defined as in (4.1)-(4.4) followed by analytic continuation into 
C\£ and the polynomial matrix pencils T n ,T,, Q i,n-i,s, 0 2 ,n- i,e, and 7Y n+ i )E are 
introduced as follows: 

(?) o) is an mxm monic matrix pencil of degree n, that is, E n ,i:{-,x o) is 

of the type 


•7>»,s(z,ao) = -^o.s^o ) = I m , z € C (4.17) 


t=o 


and 


2R 


. T n s(-,xo) is assumed to be an mxm. Herglotz matrix. (4-18) 

1/Z ’ 


'2n+l 


Hence E n< Ti{-, xo) is a self-adjoint (in fact, strongly hyperbolic) pencil, 
E n ^(z, £o)* = E n ^{z, x 0 ), z£C 

and (j/2)i?^i^„,E and 2ii?2n+i-7 r I )~i; admit the Herglotz representations 


(4.19) 


E„,n( Z , Xo) = 


2fi 2?l+1 (^)V2 


2tt 7 s 7?2„ +1 (A) 1 / 2 


•7>i,e(A, a.’o)-r——, 2 GC\S, 


A — z 


(4.20) 


*i?2n+l(2) 1/2 l C n,E(^,a:o) 1 

1 


dX R 2n+ 1(A) 1 / 2 ^„ iS (A, a-o)" 1 


71 Js 


1 


A 


A - z 1 + A 2 


A 


+ r s ,o(a:o) ^ Y,( z - Hk{ x o)) 1 r Eifc (a;o), 


(4.21) 


fc =i 


z G C\{£ U {/ifc(a:o)}i<fe<A}i 


where 

r E ,o(*o) = r E ,o(xo)* G C mxm , 0 < r Ejfe (cr 0 ) G C mxm , 1 <k<N, 

N n 

Y rank(r Eifc (a;o)) < mn, /J. k (x 0 ) G (J [E 2 j-i, E 2j ], 1 < k < N. (4.22) 

k —1 3 =1 

In fact, there are precisely to numbers ^k{xo) in [E 2 j-i, E 2 j\ for each 1 < j < n, 
counting multiplicity (they are the points z where E n ,i:(z, xo) is not invertible). 
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(ii) Given these facts we now define 


N 


Sl,„-l,s(*,*o) = (£T3 <a '“ ) 


k=1 


z - fJb k (x o) 


N 


' k=1 


and 


^n+l,E(-2’) *^o) — -^2n+l (-2')*^'n,X)(-^j 3?o) 

W £fc(x 0 ) 


jT E>fe (xo)Jl p „ i E;(^, xo), 

(4.23) 

' £fc(x o) r , A 
[ 2 - Pk(x 0 ) J 

(4.24) 

z € C\{/ifc(x 0 )}i<fe<jv, 

(4.25) 


(4.26) 


g^- Mfc (x 0 ) rE ’ fe(a:o) )- Fn ’ E( ^ a;o) ( 


AT 


£t(x 0 ) 


«=i 


w(*o) 


rs,£(a;o) ), 


£ C\{/Xfc(a;o)}i<fc<Ar. 


Lemma 4.3. Lei 2 ; € C\{^fc(:ro)}i<fc<Ar • 0p,ra-i,E(’> Xo), P = 1,2, are m x m 
polynomial matrix pencils of equal degree at most n — 1 and TL n + i,e(',xo) *s a 
self-adjoint m x m monic matrix pencil of degree n + 1. Moreover, the following 
identities hold. 


G 2 ,n-l,z(z,X 0 )* = Gl,n-lM Z ’ X o), (4.27) 

d- n ,T.{z, xo)Gi, n -i,i:(z, xo) = l/2,n-i,s(^, xo^n.E^, Xo), (4.28) 

Hn+l,E(z, So)^2,n-l,E(-2, Xo) = </l,n-l,s(z, Xo)7d„+i,S ( 2 , Xo), (4.29) 

A,e(-2,Xo)'H„ + 1 i e( 2,X 0 ) - l/2,n-l,E(^,X 0 ) 2 = i?2n+l(-s)^m, (4.30) 

W„ + 1,e(^,X 0 )1 C '„,e(2;,Xo) - 1/1,„-1,e(^,X 0 ) 2 = i?2n+l(-)^m- (4-31) 

Proof. The identities (4.28)-(4.31) are obvious from (4.23)-(4.26). Similarly, (4.27) 
is clear from (4.23)-(4.25) and (4.19). By (4.21), one infers 

lim T-z,k(xo)3 r n,n(z,Xo) = 0 = lim P n ,s(z,x 0 ) r Eife (x 0 ) (4.32) 

Z-*[lk(xo) Z—>tilc(xo) 

and hence G P , n - i,e(‘> Xo), P = 1,2, are polynomial matrix pencils of degree at most 
n-1. By (4.27),’ 

deg(^i >n _i >E (-, x 0 )) = deg(^ 2 ,n— i,e(*, x 0 )) < n - 1. (4.33) 


Next, using (4.21) again, one notes that 

ii?2n+i(Mfe(^o)) 1/2 rs.fe(xo) = -r s> fc(xo)[(d/d2;)J 7 „ iS (/afc(xo))]r s , fc (xo), (4.34) 

1 < k < N 

and thus, combining (4.26) and (4.34), 

res z=Mfe ( Xo )W n+ i iE (2;,xo) (4.35) 

= »^ 2 n+i(Mfe(xo)) 1/2 r s , fe (xo) + T- Stk (x 0 )[(d/dz)J r n) s{nk(xo))]Ts t k(x 0 ) = 0, 

1 < k < N. 

Hence, H n + i jE is indeed a polynomial matrix pencil of degree n + 1. □ 


In fact, Hn+i ,e is a strongly hyperbolic pencil as shown in Theorem 4.8. 

In the following it will be convenient to use the following set of assumptions. 



MATRIX SCHRODINGER OPERATORS WITH FINITE-BAND SPECTRA 


17 


1 / 2 

Hypothesis 4.4. Let m,n € N. Define R 2n +i as in (4.1), (4.2) and R 2n+1 as 
in (4.3), (4.4) followed by an analytic continuation to C\£, with £ introduced in 
(4.5). Moreover, let the polynomial mxm matrix pencils Q i,n-i,s(‘, £o)> 

52,h-i,s('i^o)) and 7t n +i,s( , ) xq) be defined as in (4.17), (4.18), (4.23)-(4.26). 

Next, we introduce 

M±,x(z,x 0 ) 

= ±iR 2 n+i{z) 1/2 iFn^(z, Xo)- 1 - Q\,n-\,n{z, x 0 )T n ,Y.{z, Zo) _1 (4.36a) 

= ±iR 2n+1 {z) 1/2 J 7 n ^(z, x 0 ) 1 - T n ,Y.{z, x 0 )~ 1 g 2 , n - 1 ^(z, x 0 ), (4.36b) 

2 ; G C\{£ U {Abc(zo)}i<fc<iv} 

and 

■A/±,e(z, Xo) = M- t v(z, x 0 ) ± M + ,z(z, x 0 ), z G C\{£ U {pk(x 0 )}i<k<Nj- 

(4.37) 

We also introduce the open interior £° of £ defined by £° = U"=o(Sy>i$y+i) u 
(E 2n , oo). Then one verifies the following fundamental facts. 

Theorem 4.5. Assume Hypothesis and let z G C\{£U{/ife(a;o)}i<fc<Ar}- More¬ 
over, introduce the 2 in x 2 m matrix My(-,xq) as in (4.16), (4.37) and the in x m 
matrices M±,s(-,x o) as in (4.36). Then, 

(i) ±A4±,s(-, xo) are in x in Herglotz matrices with representations 

±M±,y(z,x 0 ) = ^ J dAi?2n+i(A) 1/2 ^ r „,E(A,a;o) _1 “ 1 ^\ 2 ) 

+ rs,o(a:o)^^- k \ °! rs,fc(^o) ; (4.38) 

f!f{ Z ~ Tk\xo) 

z G C\{£ U {Mfc(^o)}i<fe<iv}- 

Moreover, M±,y(',x o) are the half-line M-matrices associated with self-adjoint 
Schrodinger operators H± Xo Y , in L 2 ([x o,±oo)) m , with a Dirichlet boundary condi¬ 
tion at the point x$ and an m x m matrix-valued potential Qy satisfying 

Qy = Qh G L 1 1 oc (R) mxm n C°°(K\{xo}) m><m , (4.39) 

given by 

Hg tX0> v = -I m -^ + Q s , (4.40) 

clom(i7| ) ia . o>s ) = {gG L 2 ((x 0 ,±oo)) m \g,g' G AC([x 0 , c]) m for all c ^ x 0 ; 

lim (/(x 0 ± e) = 0; (- g" + Qsff) G T 2 ((a; 0) ±oo)) m }. 

ej.0 

(ii) The differential expression Cy = —Dmfjfiz + Qs is in the limit point case at 
it oo. 

(Hi) The matrix A4 y(-,Xq), defined by 

My(z,x o) = (Af 5 ;, M (z,io)) 1 < M < 2 (4.41) 

_ / A^±,E(z,a;o)A7^ 1 E(z,a;o) -1 At^,E(z,a:o) H-,s(z,xo)~ 1 Af+,s(z,xo)/2 A 
V A7 +lE (z,a;o)A7-,s(z,a;o) _1 /2 A7_,E(z,a;o) _1 ) ’ 
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is a 2to x 2to Herglotz matrix admitting a representation of the type (2.43), with 
measure Qy(‘,x) given by 


dflui X,x 0 ) 


{ 1 / 

, Hn+ i,s(A,ajo) 

— ^2,n-l,s(A,lEo) \ 

J 27Tfl 2n + 1 (A) 1 /2 \ 

\o, 

v — ^i,n-i,s(A,a:o) 

^n,E(A,a;o) / 


dX, 


A e E°, 

A £ R\E. 

(4.42) 


In addition, X4 y{',xq) is the Weyl-Titchmarsh M-matrix associated with the self- 
adjoint Schrodinger operator Hy in L 2 (R.) m defined by 


d 2 

Hy = -Zmj-2 + Qs, (4.43) 

dom (Hy) = {g £ T 2 (K) m | g,g' £ AC loc (R) m ; (-g" + Q s g) £ L 2 (R) m }. 


(in) ifs /ias purely absolutely continuous spectrum E, 

spec(i? s ) = spec ac (i? s ) = E, spec p (7J Sri ) = spec sc (i?s„) = 0, (4.44) 

with spec (Hy„) of uniform spectral multiplicity 2m. 

Proof. The representations (4.38) for A4 ± i e immediately follow from combining 
(4.21), (4.23), and (4.36a). These representations also prove that ±A4± > s(-,xo) 
are m x m Herglotz matrices (cf. Theorem 2.3). Combining (4.38) with Theorem 
2.4, taking into account Lemma 8.3.2 in [65], then yields the properties stated for 
Qy- That Cy = + Qy is in the limit point case at ±oo can be proved 

in analogy to Wienholtz’s proof [109] of a result originally due to Povzner [84], 
reproduced as Theorem 35 in [41, p. 58]. The corresponding details will be presented 
in [22]. Equation (4.41) follows from (4.16), (4.23)-(4.26), and (4.36). Relation 
(4.44) follows from the explicit formula (4.42) of the spectral measure. In particular, 
the support property supp(fls) = E of the measure fis in (4.42) proves spec (Hy) = 
E, etc. The uniform maximum spectral multiplicity 2m then follows from the fact 
that rank(df!s/dA) = 2 m on the interior E° of E. □ 


At this point we cannot yet infer continuity of Qy at the boundary point xo- We 
will subsequently return to this issue in Theorem 4.8. 

In the following we will apply these facts to our concrete class of matrix-valued 
Schrodinger operators discussed in Theorem 4.5. In order to find the corresponding 
Weyl-Titchmarsh matrices X4±,y(z, x), we need some preparations. We denote by 
if±x{z, x, xq) the Weyl solutions (2.31) associated with Qy, that is, 

i>±,T.{z, x,x 0 ) =9y(z,x,x 0 ) + 4>Y(z,x,x 0 )M± t Y(z,x 0 ), z € C\E, (4.45) 


where, in obvious notation, 9y(z, x, xo), 4>y(z, x, xq) denote the fundamental system 
(2.4) corresponding to Qy- Then straightforward computations of the right-hand 
sides of (2.46)-(2.49) (taking into account (2.16), (2.17), (4.28), and (4.31)) yield 

My(z,x) = (14s, m (2,i)) 1$m$2 (4.46) 

Hn+l^iz^) -<? 2 , n - l ,y(z,x)\ y 

~Gl,n-l,v(z,x) Fn,Y.{z,x) )’ ' ’ 


2R 2n+1 ( z y/ 2 
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where we abbreviated 


n,s(z, 

X) 

II 

M 

> a;o)^ r n,E(z, x o )0e(z, x, x 0 )* 





+ </>e(z, 

, X, X 0 )Wn+l,E(z, Xo)(/>E(z, X, 

Xo)* 




- (f>Y.{z, 

, x, Xo)^l,n-l,s(z, X 0 )6»s(z, X 

,xo)* 




- 6»e(z, 

X, Xo)t/ 2 ,n-l,s(z, X 0 )(/>e(z, X 

,x 0 )*, 

(4.47) 

—1,e(z, 

X) 

= —9'^(z, 

x, xq)!F n,s(z, x 0 )6>e(z, x, x 0 ) 

* 




— 0e( 2 : 

, X, Xo)7f n-{- 1,E (z, Xo)0E (z, X, 

Xo)* 




+ ^e( 2 : 

, x, X O )01 ) „-1 ,e(z, x 0 )6»e(z, X 

,Xo)* 




+ ^s(z, 

x, Xo)02,n — 1 ,e(z, Xo)(/>E (z, X 

,Xo)*, 

(4.48) 

^/2,rt—1 ,e(z, 


= -0s(z, 

X, X 0 )^ r „,E(z, X 0 )9'- S (z, X, Xo) 

* 




- ^e(z ; 

, X, X 0 )W„ + 1,E (z, X 0 )(/>s(z, X, 

Xo)* 




+ </>e(z. 

, X, x o )0i,„to,e(z, x o )0s(z, X 

,xo)* 




+ 0e(z, 

x, X 0 )f/ 2 ,n — 1 ,E (z, X 0 )(/>s (z, X 

,Xo)*, 

(4.49) 

h~tn-\-l ,E (z, 


= #s(z,X ; 

, Xo)^ r n,s(z, X 0 )6^(z, X, Xo)* 





+ ^e( 2 : 

, X, X 0 )H„ + 1,e(z, Xo)</4(z, X, 

Xo)* 




~ 4>'t. (z : 

, X, x o )0i,„_i,e(z, x o )0s(z, X 

,Xo)* 




— ^eG 2 ) 

X, Xo)f/ 2 ,n-l,E(z, X 0 )(/>s(z, X 

,xo)*. 

(4.50) 


Considerations of this type can be found in [65, Sect. 8.2] in the special scalar case 
to = 1 and in the matrix context to G N in [101, Sect. 9.4], 

Differentiating (4.47)-(4.50) with respect to x (taking into account that 9" = 
(Qs — zX m )9 , cj>" = (Qs — z)(j>) then yields (4.51)-(4.58) below. Alternatively, 
these results directly follow from Lemma 2.9 identifying g and (i/2)R 2n ^_ 1 tF n , g p 
and (i/2)R 2 n+iG P ,n-i, = 1,2, and f) and (i/2)R 2 ^ 1 'H n +\, respectively. 

Lemma 4.6. Assume Hypothesis 4-4 an d iei (z,x) G C x R. Then 


F'n,Y, — — (£l,n-l,E + G'l.n I.e)- (4-51) 

^l,ra-l,E = '(Qs — zIm)tFn,S ~ H n + 1,E (4-52) 

= {-K,n + ^h.eQe - Qe^,e)/2, (4.53) 

^l,n- 4 ,E = (Qe - zl m - Qs T n ,Y. + 5i,t»-1,e2e — QE01,ra-l,E, (4.54) 
02,ra-l,E = ti,e(Qe - Zl m ) — W n +1,E (4.55) 

= s + Qs^„,e - ^„,E Qs)/2, (4.56) 

^2,n-l,E = — 2 J^ )S (Qe — zX TO ) - + Qs&,n-1,E — &,i»-1,eQe, (4-57) 

^n+l,E = '5l,n-l,E(Qs - zX m ) — (Qs — zJ m )t/ 2 ,n-l,E, (4.58) 

W n+1 ,E = [^",E - ^,e(Qe - Zl m ) - (Qs - zX m )^„,s]/2. (4.59) 


In particular, one also verifies the following facts from (2.56)-(2.63). 
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Lemma 4.7. Assume Hypothesis 4-4 an d iei ( z ,x) G C x R. Then 

4Fn, e(z, x)* = E n ,Y.{ z , x), H n +l,n(z, x)* = Hn+l,Tl(z, x), 

02,n-l .eCz,®)* = Gl,n-l&{z,x), (4.60) 

Fn,T.{z,x)Gl,n-l,T.{z,x) = G 2 ,n-l,T.{z,x)T n ^{z,x), (4-61) 

'Hn+l,-z{z,x)Gl,n-l,T.{z,x) = Gl,n-l,v(z,x)Hn+l,T:(z,x), (4-62) 

^n+l,s(z, x)J~ n,x;(z, x) Gl,n— 1,X3(-Z, x) /?2n+l(^)-^rrn (4.63) 

E n ,T.(Z, x)H n +l,s(z, x) - <? 2 ,n-l,s(z, X ) 2 = (4.64) 

Proof. (4.60) is clear from (4.27), the fact that E n ,T,{'i x o) and Tt n+ i,s(‘>£o) are 
self-adjoint m x m matrix pencils, (4.48), and (4.49). Similarly, (4.61)-(4.64) follow 
from elementary (but somewhat tedious) calculations directly from (4.47)-(4.50), 
invoking (2.10)-(2.17) and (4.28)-(4.31) repeatedly. □ 

Combining (2.50)-(2.54) and (4.46) then yields 
M- ;>: (z.x) 

= ±iH 2 n+i(z) 1/2 d 7 n,s(z,x)~ 1 - Gi, n -i,s(z, x)E n ,s(z, x) _1 (4.65a) 

= ±iR 2n+ i(z) 1/2 E n ^(z, a;) -1 - E n ,Y z , x)~ l G 2 , n -i,T.{z, x), (4.65b) 


One observes that for each x £ R, A4+,e(-,x) is the analytic continuation of 
M-,s(-,x) through the set E, and vice versa, 

lim_A4 +,e(A + ie, x ) = limA/f-^A — ie, x), (4.66) 

ej,0 eJ,0 

n — 1 

A G [J {E 2 j, E 2 j + \) U (E 2n , oo), If 1. 

3=0 

In other words, for each ifl, A4 +,e( - , x) and _A4 _,e(-, x) are the two branches of 
an analytic matrix-valued function AAs(-,x) on the two-sheeted Riemann surface of 
i? 2 n+i- Thus, the corresponding potential Qs is reflectionless in the sense discussed 
in [23], [40], and [56]. 

Thus, one obtains the following results. 


Theorem 4.8. Assume Hypothesis 4-4 an d let z £ C\R and ifl. Then 
(i) M±.t,{z, •) in (4.65) satisfy the matrix-valued Riccati-type equation 

M± -S (z, a;) + M± t s{z, x) 2 = Q^(x) - zX m , i f I, z f C\l. (4-67) 


Moreover, -M± ; e(z, x) in (4.65) are the m x m Weyl- Titchmarsh matrices associ¬ 
ated with H± x E on the half-lines [a;,±oo) and thus for each x £ R, A4s(z,x) in 
(4.46) is a 2to x 2m Weyl-Titchmarsh matrix associated with H-% on R. In partic¬ 
ular, AAs(-,x o) is a 2m x 2 m Herglotz matrix of Hy. admitting a representation of 
the type (2.43), with measure f Y(-,x) given by 


dfls(A, x) 


1 ( Hn+l.sfA,®) — 02 ,n-l,s(A,a;) \ n 

27ri?2n+l(A) 1 / 2 ^n,s(A,*) ) ^ 

0 , 


A G E°, 

A G R\E. 


(4.68) 


(ii) T n ,e(-,x) and 'H n +i,s(-, x) are strongly hyperbolic (and hence self-adjoint ) to x 
to monic matrix pencils of degree n and n + 1, respectively, and Gp, n - i,s(-,a;), 
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p = 1,2, are mxm matrix pencils of degree n — 1. 

(in) Q s G C°°(R) mxm . 

Proof. (4.67) is clear from Lemma 4.6 and (4.65). Since Qs G C loo (R) mxm , the 
initial value problems 

M'±(z, x) + M±(z, x) 2 = Qt.(x) — zl m , x G R, z G C\R, (4.69a) 

M±(z,x 0 ) = M±^(z,x 0 ), (4.69b) 

with M±.y.(z,Xq) given by (4.36), has a unique solution. Since at x = Xo this 
solution coincides with the Weyl-Titchmarsh M -matrix M±x(z, x)\ x = Xo in (4.65) 
(using the initial condition (2.5) in (4.47)-(4.50)), M±^(z, x) represents the Weyl 
Titchmarsh matrices associated with on the half-lines [a;,-too). This proves 
part (i). By the known leading asymptotic behavior (2.33) of A4±,s(-,x) (valid 
for each xo G R, see [21]) and that of the diagonal Green’s matrix Gz(-,x,x) = 
A4e, 2,2(')X) of Hy. as | z | —> oo, J~n,Y .(*, x) and 'H n +i,T .(•, x) are monic matrix pencils 
of degree n and n + 1, respectively, and G P , n - i,e(-, x), p = 1,2, are m x in matrix 
pencils of degree n — 1. Since the diagonal blocks of each Herglotz matrix are 
also Herglotz matrices, one concludes that (i/2)Rffn+\J : n,T; and (i/2)R~ 2 n+l^n+l.S 
are Herglotz matrices. By Corollary 4.2 this then proves that T n p, and H n + i,s 
are strongly hyperbolic pencils and hence item (ii) holds. As in Theorem 4.5, 
Qs G C°°((— oo, x) U (x,oo)) mxm . Since x G M is arbitrary, this proves (iii). □ 

It should be emphasized that the construction of Qs along the lines of Section 4 
in the scalar case m = 1 is due to Levitan [62] (see also [63], [64], [65, Ch. 8], [67]). 

5. Trace Formulas and Connections with the 
Stationary Matrix KdV Hierarchy 

In this section we introduce the stationary matrix Korteweg-de Vries (KdV) 
hierarchy and show that the class of finite-banc! potentials Qs constructed in Sec¬ 
tion 4 satisfies some (and hence infinitely many) equations of the stationary KdV 
equations. We also introduce trace formulas for KdV invariants. 

In order to extend the recursive approach constructing KdV Lax pairs in the 
scalar (Abelian) context to the present matrix-valued (non-Abelian) setting, we 
focus on an efficient approach introduced by Dubrovin [28] (in the scalar case m = 
1). Recalling (2.50)-(2.63) and Lemma 2.9, we state the following matrix-version 
of Dubrovin’s generating function approach to higher-order Lax pairs. 

We start by introducing the following hypothesis. 

Hypothesis 5.1. Fix m G N, suppose Q = Q* G C'°°(R) mxm and introduce the 
differential expression 

d 2 

C = —X TO + Q, x G R. (5.1) 

Suppose C is in the limit point case at ±oo and introduce the corresponding self- 
adjoint operator H in L 2 ( R) m by 

d 2 

H ~ - lm di 2 + Q ’ 

dom(H) = {g G L 2 ( R) m \g,g f G AC loc (R) m ; (-g" + Qg ) G L 2 ( R) m }. 


(5.2) 
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Given Hypothesis 5.1, we introduce the generating operator P z by 

P z = ( 0 ( 2 , -)f- + 02 ( z , •))(£ - zX m )~ l , z G C\R, (5.3) 

dom(P z ) = [J ker (C-El m ), 

-EeC\{RU{z}} 

where ker(£ — EX m ) denotes the algebraic nullspace of C — EX rn (as opposed to 
the functional analytic nullspace in L 2 (R) m ) and (£ — zX m ) _1 acts in the obvious 
manner by 


(£ - zl m ) 1 tp = (E - z ) 1 tp, tp G ker(£ - EX m ). (5.4) 

The precise operator theoretic properties of P z will be irrelevant in the following. 
Lemma 5.2. Assume Hypothesis 5.1 and let z G C\R, iGl. Then, 
d 


(flObO^;+fl 2 {z,-)){C- zXm) \ 
= - 


C 


UEeC\{RU{z}} ^ er (^ EXm) 


UfieC\{RU{2}} EXrn) 

Proof. Fix ^ G ker(£ — EX m ) for some E G C\{M U {z}}. Then one computes 
d 


(5.5) 


+Q2(z, ■)){£-zXm) \£p 
= (E-z)- 1 {[2q ' 2 + g"+sQ-Q&W 
[02 + 2g ( Q — zX m ) + qQ + 02 Q — Q 02 + 2 (z — E)g ]ip} 

= -2flV, (5-6) 

using (2.68) and (2.72). □ 


A variant of Dubrovin’s idea of a generating operator for KdV Lax pairs was 
also used by Olmedilla, Martinez Alonso, and Guil [79]. Their approach, however, 
focuses on formal operator expansions and formal pseudo-differential operators. For 


a different approach we refer to [27, Ch. 15], [32], 

Next, we recall that 

g{z,x) = [M-(z,x) - M + (z,x)]~ l , (5.7) 

01 (• z,x) = \[M-{z,x) + M+{z,x)][M-{z,x) - M+(z, a;)] -1 , (5.8) 

02 (z,x) = ^\M-(z,x) - M+iz^^^M-iz^) + M + (z,x)\, (5.9) 

f ){z,x) = M±{z,x)[M-{z,x) — M+{z, x )] -1 M^{z, x), (5.10) 


and note that by Theorem 2.7 the right-hand sides of (5.7) and (5.9) admit as¬ 
ymptotic expansions in cones avoiding the spectrum of H. In particular, one thus 
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obtains the asymptotic expansions 


OO 


g{z,x) 

\z\ —KX) 

-±^J2^( x )z~ k , Wo (x)=l m , 

(5.11) 


zec e 

k—0 


9p(z,x) 

\z\—>oo 

27I/2 ^ v &p,k(%)z ? e Pl o(*) — ^5 p — 2, 

(5.12) 


zdC e 

k—0 


t)(z,x) 

\z\— »oo 

1/2 °° _ 

(5.13) 


zec e 

k—0 



for some coeficients fftfc, & p ,k, P = 1,2, and 9) k , which are universal differential 
polynomials in Q by Remark 2.8 (i). Explicitly, one obtains 


&o=X m , ^r = IQ, m 2 = -iQ" + |Q 2 , 

(5.14) 

©i,o = -iQ', © 1,1 = ^2'" - |(2 2 )' - f Q'Q, 

(5.15) 

© 2,0 = -iQ', © 2,1 = tqQ!" - |(Q 2 )' - IQQ ', 

(5.16) 

$o=2m, $i = -§Q, ^2 = |Q"-|Q 2 , 

(5.17) 

etc. 



Motivated by Lemma 5.2, we now introduce the in x m matrix-valued differential 
expressions P 2 k+i by 

fc , 

p2k+i = '52(*t-r + 8*’t) ck - e ’ (5-is) 

£—0 x 

In analogy to Lemma 5.2 one then obtains the following result. 

Lemma 5.3. Assume Hypothesis 5.1 and let z G C\R, i£l. Then, 

[V 2 k+i,C] = 2 m' k+1 , k G No. (5.19) 

Proof. Assuming ip G ker(£ — zl m ) one computes 
k 

[V 2k+1 , c] Ip = Y, z k ~ £ ( 2 &' 2/ + % + Q - Q&e) P’ 

e=o 

k 

+ J2 z k ~ e {& 2 ,e + 2 %(Q - zl m ) + ft e Q' + & 2 ,eQ - 0 & 2 ,i)P 
£—0 
k 

= z k ~ t (2%' 2/ + <R" +<R e Q - QSktW - 

k 

+ J2 z k ~ e {&ie - 2 ^+i + 2^Q + %Q! + © 2 , e Q - Q© 2,i)p + 29t' fc+1 V> 

i =0 

= ^' k+ iP- (5-20) 

Here we used 91 q = 0 and 

2©^ + m'/ + m e Q- Q&e = 0, l G N 0 , 

02,/ - 2^+i + 2 *k t Q + m e Q' + - Q® 2 / = 0, £ G N 0 , 


(5.21) 

(5.22) 
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which follow from inserting the asymptotic expansion (5.11) and (5.12) into (2.68) 
and (2.72) (which is permitted by Theorem 2.7). Relation (5.20) implies (5.19) 
since V-ik+i and C are m x in matrix-valued differential expressions of finite-order 
while LLec ker(£ — zl m ) is an infinite-dimensional space of C°°(M) m -functions. □ 


Introducing 


k 

~P2k +1 = y^Cfc-iT^i+i) k € No, (5.23) 

e=o 

where 

C C, co = 1 (5-24) 

denotes a set of constants, the pairs {fP 2 k+\i £)> & G No, by definition, represent the 
Lax pairs of the (matrix-valued) KdV hierarchy. More precisely, varying {c^l^eN C 
C, the set of evolution equations, 

j t C-[V2k+ 1 ,£]= 0 , fceNo, (5.25) 

or equivalently, the set of equations, 

k 

KdV fc (Q) = Q t -2^c fc _^ +1 (Q,...) =0, fceNo, (5.26) 

^=o 


represents the (matrix-valued) KdV hierarchy of evolution equations. The corre¬ 
sponding stationary KdV hierarchy, characterized by 


Qt = 0, or equivalently, by [P 2k + i,X] = 0, fee N 0 , (5-27) 


is then given by 


s-KdV fc (Q) 


-2^c fc _^ +1 (Q,-...) =0, he N 0 . 


«=o 


(5.28) 


Remark 5.4. By Remark 2.8, each 8SKt is a universal polynomial in Q and its x- 
derivatives and occasionally we slightly abuse notation and indicate this by writing 
9^(Q,...) for y\e(x), 93( +1 (<2,...) for fft( +1 (a;), etc. Explicit formulas for the 
differential polynomials 91 e(Q ,...) were derived for instance, in [8] and [83]. 
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In order to make the connection with the finite-band formalism of Section 4 we 
now recall 


2iW(z) 1 /2 


Fn,x(z,x) = [A4_ e(z, x) - M + t -z{z,x)\ i , 


(5.29) 


2I? 2 „ +1 ( 2 )i/2 


01,n-lM Z ,x) 


= -[M-,z(z,x) + M+,-z{z,x)][M-&{z,x) - M+^(z,x)\ \ 

G2,n-1,T.{Z,X) 


2R 2n+1 (z) 1 / 2 


= -\M-^{z,x) - M+v(z,x)\ l [M-,t,(z,x) + M+,-z{z,x)\, 
r^n+l,E (z,x) 


2R 2n+1 (z) 1 / 2 ' 

= M± i e(z,x)[M- i e(z,x) - M+x{z,x)Y 1 M t ,z(z,x), 


(5.30) 


(5.31) 


(5.32) 


and note that (5.29)-(5.32) admit expansions convergent in a neighborhood of in¬ 
finity. In particular, 


1 1 00 ^ 


kh<» z 1 
zec e k ~° 


1 i 00 _ 

, , ) i/ 2 Sp,n-l,s(z,x) = 7372 ^ ®p,kM x ) Z k i 


R 2 n+i{z) 1/2 


PHoo Z 1 / 2 
zec e k 0 


93o.s(x)=J m , (5.33) 

®p,o,E(a:) =0, p= 1,2, 
(5.34) 


iW(*) 1/2 


oo 

= -777 V'£)fc,E(:r) 2 -fe , fio,v{x)=l„ 

1 2 1 —>00 Z L ' Z 


(5.35) 


for \z\ sufficiently large. Here the coefficients and © 2 ,fe,s are the universal 
differential polynomials 94*, = 94fc(Qs> • • •) and 0 2 ,fc = ® 2 ,fe(Qsi • • •) hr (5.11) and 
(5.12) (with Q replaced by Qs)- We also recall 


^nM Z i X ) = (x)z e , JZq^x) =l m , (5.36) 

t=0 
n —1 

Gp,n-l,s{z,x) = ^ Gp,n-l-t,Y.{x)/, p= 1 , 2 , (5.37) 

t=0 

n 

H,i+i^(z,x) = '^2'H n+ i-e^(x)z e , Ho,n(x) = l m . (5.38) 

e=o 

Since we seek the connection between the set of coefficients 94fc,s, ® P ,fc,s, and 
•Ttj.Ej Gp,k,Y.i 7~tk ,Ej we next consider the following elementary expansions. Let 

77 G C such that |?j| < min{|.Eo| _1 , • • • , l-E^nl -1 }- (5.39) 



26 


F. GESZTESY AND L. A. SAKHNOVICH 


2 n \ —1/2 oo 

II (i ~ EfJI ) ) = d k{E)ii k , 

e=o ' fc=o 


(5.40) 


where 


co (E) = l, 


«(£) = E *«"• 


JO t ■ ■ ;j2n—0 

ioH-hj2n=fc 


2 2fc Oo0 2 --'(j2n!) 


(5.41) 


The first few coefficients explicitly read 


cq(E) — 1) c 1 (E) = -Y,E e , 


02 (E) — - ^ Ej 1 Ee 2 + ^ y/ E 2 , etc. 
* 1 / 2=0 0 *=o 

*i<* 2 


Similarly, one has 


2n \ 1/2 00 

II (1 - E w )) =^2 ck(E)ii k , 

*=o ' k -0 


(5.43) 


where 


co(E) = 1, 


( F \ _ _ _ (2jo)! • • • (2j 2n )! _ j 0 j 2n 

;Uj y _g2 2fc (j0 ! ) 2 • • • (j2n!) 2 (2jo — 1) ' ' ' (2j2n — 1) °"‘ 2 "’ 

J 0 v J 2 n—U 
ioH-bj 2 n=fc 


fcel (5.44) 


The first few coefficients explicitly are given by 


co(E) = 1, 0!^) = --^^, 


^ 2n ^ 2n 

2 (E) = T ^2 E (l E #2~ Q^2 E e> etc - 


«1,<2=0 

^1<^2 


(5.45) 
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Lemma 5.5. Assume Hypothesis 5.1 and let ieR. Then, 


t 


Ei,y.(x) = Y j ce-k{E_)iRk,z(x), 
fc=0 

i 

l = 0,. 

. .,n, 

(5.46) 

^A.e,T.{x) = Y. ce-k(E_)Ek,Y.(x) , 
fc=0 

t 

t = 0,. 

..,n, 

(5.47) 

Qp,t,T.{x) = Y< Ct-k{.E_)®p,k+l,Y.{x 
k =0 

t 

), e = 

0,..., n — 1, p= 1,2, 

(5.48) 

®p,<+i,e(*) = Y', ce-k {Ef)Qp,k,Y. (x ), 

fc=0 

e 

£ = 0, 

..., n — 1, p= 1,2, 

(5.49) 

He, b(x) = Y c e-k(E)£) k ,x(x), 
k=0 

e 

e = o,. 

..,n+l, 

(5.50) 

fieMx) = Y^-k(R)HkMx), 

k =0 

e = o,. 

• • ,n + 1. 

(5.51) 


Proof. It suffices to prove (5.46) and (5.47) and so we omit the analogous proofs of 
(5.48)-(5.51). Since for \z\ sufficiently large, 

n 

Z~ n F n ^(z,x) = (x)z~ f 

£=0 


= z- n ~^R 2n+1 (*) 1/2 £ (x)z~ t 
r=o 

OO OO 

= Y, c k(E)z- k J2^( x ) z ~ e 


k —0 


1=0 


-l 


= X] ( y c (.-k{Ef}‘5\k,Y.{x) 

1=0 ' k =0 

and hence (5.46). Equation (5.47) is then clear from (5.46) and 

k 

Y Ck-e(E)ce(E) = 6k,o, k e N 0 . 

e=o 

The latter follows from multiplying (5.40) and (5.43), comparing coefficients of 

77 fc . " □ 


(5.52) 


(5.53) 


Given these preliminaries we can now state the principal result of this section. 

Theorem 5.6. The self-adjoint finite-band potential Qs € C°°(R) mxm , discussed 
in Theorems f.5 and j.8, is a stationary KdV solution satisfying 

n 

s-KdV„(Q s ) = -2]Tc„_r(E)^ +1 (g E ,...) = 0, (5.54) 

t=o 

with ce(E_) defined in (5.44) and Tie+i the universal differential polynomials (with 
respect to Q) in (5.11). 
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Proof. Introducing the m x m matrix-valued differential expression V 2 k+i ,s by 


Vin+l, E = Y. C n -l(E)V‘M + i = ^ {fn-t, e (')^ + f? 2 ,ra-l-^E(’)) ^ 
£=0 £=0 

(cf. (5.46), (5.48), and (5.23)), one computes for ip £ ker(£ — zX m ), z £ C, 


(5.55) 


V2n+l,Y.i> = Tn,Y.(z, x)lf' + <?2,71-1,e(<2, z)V’> (5.56) 

and hence, 

PWl.E, C]fl> = [2&,„-i,e + ^"s-^.sQs - Qe^.e]/ (5-57) 

+ [£/2,n-l,E + 2^ r n,s(2s - zX TO ) +^7 i,eQe +f/2,n-l,sQs — Qs&.n-l.s]^ = 0 

by (4.56) and (4.57). Since z £ C is arbitrary, this implies 


P^n+l.E) £] = 0, 

completing the proof by (5.27), (5.28). 


(5.58) 

□ 


Next, we turn to a discussion of trace formulas for the finite-band potential Qe 
in terms of (matrix) roots of T n ,e and W n +i,E- 

Theorem 5.7. Let (z,x) gCxl and assume Qe to be the self-adjoint finite-band 
potential discussed in Theorems f.5 and f.8. In addition, let the monic self-adjoint 
matrix pencils andTt n + i,e be given by (4.47) and (4.50). Then x) and 

TL n +\,n{-,x) are strongly hyperbolic pencils and hence admit the factorizations 


Xn,s(z, x) — ( zT m lAnixfj^zTjyi lA ri —\{xf) t ’ t {zT ra - 

-Ui(x)), 

(5.59) 

H n + l,n(z,x) = ( zl m - V n (x))(zlm - V n -l(x)) ■ ■ • ( zl m - 

- Vo (a:)), 

(5.60) 

with 



spec(V 0 (a;)) = spec(W„+i,E) n A 0 (H„ + i,s) C (-00 ,E 0 ], 


(5.61) 

spec (Uj(x)) = spec(^„,E) n A^X^e) C [E 2 j-i,E 2 j], 1 < 

j < n, 

(5.62) 

spec(Vj(a:)) = spec(W n +i,E) n A 0 (H„ + i,e) C [E 2j -i,E 2j \, 

1 < j < n. 

(5.63) 

Moreover, one obtains the sequence of trace formulas 



n k 

Y Uj k (x) ■ • • u n (x)U h (x) = Y Ck-i{E)&eMx), 

1 < k < n, 

(5.64) 

jl,j2,-.,jk—0 £=0 

ji<h<-<jk 

l. 



K, 

0 = Y. c k-e(K)^e, s(x), 

e=o 

k > n + 1, 

(5.65) 

n k 

Y v Jk(x) ■ • • Vn i x Wjo( x ) = y^c fc -^(£)^,E(x), 

1 < k < n, 

(5.66) 

jo,ji,..,jk=o e=o 

jo<ii<-<3k 



k 

0 = Y c k-e(E)fieMx), 

k > n + 1. 

(5.67) 


£=0 
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In particular, in the special case k = 1, Qz satisfies the trace formulas 

2 n 


x ^=0 ' j=l 

✓ 2n \ n 

= -( Y J Et)x m + 2Y J Mx). 


e=o 


k -o 


/n addition, one obtains 


( r ) ^ ^oo/p\mxm pi j^oofnp^mxm 


GiV G G 


for all r G Nq. 


(5.68) 

(5.69) 

(5.70) 


Proof. By (5.29) and (5.32), i/R 2n+ i J>»,e(-,x) and V-^n+i^+bsO, a: ) are Her- 
glotz matrices and hence (4.14) and (4.15) apply. In particular, tF n ,T.{-,x) and 
Tdn+ i,s(-,a;) are strongly hyperbolic pencils. Since both are monic, the factoriza¬ 
tions (5.59) and (5.60), as well as (5.61)-(5.63), hold by Theorem 3.5 (ii). By (5.33) 
one infers 

n / n \ 

z~ n T n ,Y. {z, x) = U jk (*)••• U j2 (x)U h (x) J z~ k 

fe=o ' h,h,—,jk=o 
ji<h<—<jk 


= z n {1/2 ' , R 2 n+i(z) 1/2 V'^fc,s(a;)2 
l*|-*°0 fGn 

zec £ k ~° 

/ 2 n \ 1/2 oo 

= na-w-)) 

*€C. *=° 7 fe=0 


-fc 


—k 


oo / k 


zec £ 


= ( y^ci(£)*Kfc-i,s(s) 

p| ^°° 


/ c — 0 x ^=0 


^ — k 


(5.71) 


Comparing coefficients 2 ; fc , fc G N, then yields (5.64) and (5.65). In the special 
case k = 1 one infers 


n .. / 2 n \ 

- £ U, (x) = -- ( £ ^ J T m + &i, E (a-) 


(5.72) 


j=i x e=o 

and since = <2 E /2 by (5.14), the trace formula (5.68) for Q E results. (5.66), 
(5.67), and (5.69) are proved analogously. By (5.62), U) G L°°(IR) mxm , 1 < j < n 
and hence 

91fc G C ,00 (M) mxm D L°°(IR) mxm , fc G N 0 (5.73) 


(since Q E G C 00 (M.) mxm by Theorem 4.8 (iii)). An analysis of the recursion relation 
(2.49) for M± t k combined with (5.7), (5.11) then proves that 9T/ C is of the form 

&=4Q (2fe “ 2) +^ fc (Q (2fe - 4) ,...), fc> 2, (5.74) 

with dk G K. appropriate constants and abbreviating a differential polynomial 
in Q which contains Q( 2fc_4 ) as the highest derivative of Q. Hence one infers 
(5.70). □ 
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The factorizations (5.59), (5.60), eigenvalue distributions (5.61)—(5.63), and trace 
formulas (5.64)-(5.69) are extensions of well-known formulas in the scalar case 
m = 1 (see, e.g., [28], [31], [36], [38], [50], [67]). 

Finally, the property Qe € C' 0O (R) rraxm in Theorem 4.8 (iii) can be improved 
upon by using the system (4.51), (4.52), (4.55), and (4.58). In fact, writing 


n 


^n,E (z,x) = ^2,F n -t,Y. (x)z e , tFo,x(x) = l m , 

(5.75) 

e=o 


71 — 1 

Gp,n- 1,e4, X) = ^ Gp,n-l-eMx)z e , p= 1,2, 

(5.76) 

t=0 


n+1 


'Hn+lMz, X) = ^2 ^n+l-t,E {x)z £ , Hq,s{x) = Z m , 

(5.77) 




one obtains the following result. 

Lemma 5.8. Assume Hypothesis 4-4 and let (z, x) G C x R. Then the coefficients 
in (5.75)-(5.77) satisfy the autonomous nonlinear first-order system 

= — (£?m-i,e + Q 2 ,e-i,s), 1 < i < n, (5.78) 

<?i/,e = — - 7G,e).?y+i,e + ^+ 2 ,e - ?4+2,e, 0 < t < n — 1, (5.79) 

= “•^+1,e('?’i,E _ ^1,e) + ^<+2,E “ ^ + 2,S, 0 < t < U — 1, (5.80) 

?4,£ = GlJ 1.>: + £/2,*-l,E — £/1,£-2 ,e(^ 1,E — Hi,s) — (^ljE — Wl,s)!?2^-2,E j. 

l<£<n + l, (5.81) 

J>1+1,E = 0, <5p,n,E = (?p,-l,E = 0, = 1,2. (5.82) 

Moreover, tFe^, 0 < l < n, Gp,t,n, 0 < l < n — 1, p = 1,2, TLt^, 0 < l < n + 1, 

and hence Qs, are all analytic in an open neighborhood containing the real axis. 

Proof. Inserting (5.75)-(5.77) into (4.51), (4.52), (4.55), and (4.58) yields 

F’h y. = ~{Gi,e-i,s + Gz,e-i, s), 1 < t < n, (5.83) 

Sm.e = -Qe^+i,e + ^+ 2 ,E - %+ 2 ,e, 0 < < n - 1, (5.84) 

^2/,E = “-4+1,eQe + ^Y+2,E - 7^+2,E, 0 < t < U - 1, (5.85) 

74,E = + !?2,f-l,£ — (4,<-2 ,eSe + Se&,<-2,e)) 1 < £ < n + 1, 

(5.86) 

•T 7 n+l,E = 0 , 0p,n,£ = f?p,-l,E = 0, p = 1,2, (5.87) 

Qe = ^1,S - Wr,E. (5.88) 

Insertion of (5.88) into (5.84)-(5.86) yields the autonomous nonlinear first-order 
system (5.78)-(5.82). Given the initial conditions 

^,s(x 0 ), 1 < i < n, 

G P ,e, s(zo), P=l,2, 0 < £ < n - 1, (5.89) 

74,e0eo), 1 < ^ < n + 1, 

determined by (5.75)-(5.77) and the half-line Weyl-Titchmarsh matrices Af±(z, Xo) 
in (4.36), the maximal interval of existence of the solution of the autonomous system 
(5.78)-(5.82), (5.89) is all of [0, oo) and (—oo, 0], and hence all of R. (cf. [42, p. 18]), 
applying (5.70). Thus, one recovers the C ,00 (M) mxm -property Ti, G P ,e, V = 1,2, 
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He, and hence that of Qs- Moreover, since the Picard iterations are convergent 
in sufficiently small circles in C centered around each x G R (cf. [44, Sect. 2.3]), 
the unique solution obtained by these Picard iterations is analytic in each of these 
circles. □ 

We conclude with a remark that puts the construction of Qe in Section 4 into 
proper perspective. 

Remark 5.9. The simplest examples of potentials Qs described in Theorem f.8 
are of the type 

Qt.{x) = diag(gi )E (a;), •.. ,q m p{x)), (5.90) 

where qj,s, 1 < j < m, are isospectral algebro-geometric finite-band potentials 
associated with scalar Schrodinger operators in L 2 (R) with spectrum E. The next 
simplest and closely related set of such examples for Qs then will be of the type 

Qs (x) = Gdiag(gi iE (x),...,gr m>s (x))G~ 1 , (5.91) 

where U G C mxm is a unitary m x m matrix independent of x el. At first sight 
one might think that perhaps all matrix potentials Qs are of the form (5.91). That 
this is certainly not the case will be argued next. Indeed, assuming that (5.91) holds 
for some unitary m x m matrix U independent of x, one infers that 


Qs\ x ) = U diag(g££(x),..., q^(x))U 1 for all r € N 0 . (5.92) 

Since by Remark 2.8 (i) and Lemma 5.5 all coefficients of T n ^, Q p , n - i,s, P — 1,2, 
and H n +i,Y. ore differential polynomials with respect to Qs, (5.92) implies 

UT n ^.{z,x)U~ x = diag{F nj x,i(z,x),..., F n ^,m(z,x)), (5.93) 

UQ p , n -\,Y.{z, x )U~ X = diag(G p , n _i, s,i{z, x),..., G Ptn -i t s, m (z, a;)), p= 1,2, 

(5.94) 

UH n +\,s.{z,x)U~ x = diag(R n+ i ! E,i(2 ; , x),.. . ,H n+1 ^ m (z,x)), (5.95) 

UM±.s{z,x)lA ~ 1 = diag(m±,s,i(^, x),... ,m±,s,m(z,x)). (5.96) 

Consequently, one obtains for all z, z' € C+, x, x' € R, 

[Q^ ) (*),Q^' ) (a: , )]=0, r G N 0 , (5.97) 

[F n ,s{z,x),F n ,s{z',x')} = 0, (5.98) 

[Gp,n-ix{z,x),g Pt n- 1 x(z , ,x')] = 0, p= 1,2, (5.99) 

[H n+1 ,s(z,x),H n+1 ,s(z',x')\ = 0, (5.100) 

[M±.s{z, x), M± t v(z', a:')] = 0. (5.101) 

In particular, 

[IFn,Y.{z,x 0 ),Fnx{z ', x 0 )] = 0 for all z , z' G C+. (5.102) 

Thus, whenever 

[Fn,Y.{z, x 0 ), F n ,Y.{z', x 0 )] ^ 0 for some z, z' G C + (5.103) 


(which can easily be arranged for n > 2), (5.91) cannot hold for a unitary m x to 
matrix U independent of x. 
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Additional results, including extensions of Borg’s and Hochstadt’s theorems in 
the special cases n — 0,1, respectively, will appear in [9]. 
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